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FACTORIZATION INVARIANTS IN NUMERICAL MONOIDS 


CHRISTOPHER O’NEILL AND ROBERTO PELAYO 


Abstract. Nonunique factorization in commutative monoids is often studied using 
factorization invariants, which assign to each monoid element a quantity determined 
by the factorization structure. For numerical monoids (co-finite, additive submonoids 
of N), several factorization invariants have received much attention in the recent 
literature. In this survey article, we give an overview of the length set, elasticity, 
delta set, w-primality, and catenary degree invariants in the setting of numerical 
monoids. For each invariant, we present current major results in the literature and 
identify the primary open questions that remain. 


1. Motivating numerical monoids 

Factorization theory has its genesis in studying the decomposition of natural numbers 
into its prime (irreducible) divisors. In this multiplicative monoid, the Fundamental 
Theorem of Arithmetic guarantees that such a decomposition is unique. Many other 
monoids, however, have elements that fail to admit a unique factorization into irre- 
ducibles; in such cases, algebraists are frequently interested in measuring, in various 
ways, how different the plural factorizations of an element can be. 

One does not have to look far for examples. In fact, certain sub monoids of the natural 
numbers have a rich non-unique factorization theory. Numerical monoids (co-£nite, 
additive sub monoids of the natural numbers) have garnered much recent interest, as 
their factorization theories are complicated enough to be interesting while controlled 
enough to be tractable. The goal of this paper is to introduce various perspectives and 
tools for studying factorizations in numerical monoids, and to present open problems 
that are of primary interest in this held. 

Several natural examples of numerical monoids are found in the literature. Histori¬ 
cally, Frobenius’ coin-exchange problem asked what monetary values one could make 
with relatively prime coin denominations. A more endearing manifestation of numerical 
monoids is the McNugget monoid, which enumerates the number of chicken McNuggets 
one could buy by using the standard package sizes of 6, 9, and 20. As with the coin¬ 
problem and the McNugget monoid, all numerical monoids can be described using a 
hnite set of generators. We provide the formal dehnition below. In what follows, we 
let N denote the set of non-negative integers. 
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Definition 1.1. Let rii, n2,..., be a collection of relatively prime positive integers. 
The numerical monoid S generated by these integers is the set of all non-negative, 
integral linear combinations of these numbers. In other words, 

S = (ni, 77-2,..., Uk) = [aiUi + 0277-2 H-h OkUk : Oi e N}. 

Remark 1.2. This set is an additive monoid with identity 0 , and the condition that 
gcd( 77 i, 772 ,..., ^k) = 1 guarantees that there exist only finitely many natural numbers 
that are not elements of S. In fact, every co-finite additive submonoid of N can be 
generated by a finite set of relatively prime natural numbers; thus, numerical monoids 
are precisely the co-finite, additive submonoids of N. 

While generating sets of a fixed numerical monoid are non-unique, there always exists 
a unique minimal generating set (with respect to set-theoretic containment). As these 
minimal generators cannot be represented in terms of the other generators, they are 
precisely the irreducible elements of the numerical monoid. For the remainder of the 
paper, we will assume that {771,772,... , 77 ^} constitutes the unique minimal generating 
set for a numerical monoid S and that 77 i < 772 < ... < 77 ^. 

Both Frobenius’ coin problem and the McNugget problem ask which natural numbers 
arise as non-negative, integral linear combinations of their respective generators. As 
such, understanding the largest natural number that is not an element of S (named 
for the aforementioned coin problem) is often useful. 

Definition 1.3. For a numerical monoid S, the values in N \ S' are called the gaps of 
S', and F{S) = max(N \ S') is called the Frobenius number of S'. 

Example 1.4. For two relatively prime positive integers 771,772, the numerical monoid 
S' = (771,772) has F{S) = 77 i 772 — (771 -|- 772) as its Frobenius number [ 50 ]. Furthermore, 
all but finitely many elements of S' can be written in multiple ways as non-negative 
integral combinations of 77 i and 772. The simplest example is the element 771772 G S', 
which can be seen as 77 i copies of 772 or as 772 copies of 77 i. These two representations of 
77 i 772 as an element of S' are called factorizations] see Section]^ for precise definitions 
and notation. 

The McNugget monoid M = ( 6 , 9 , 20 ) has Frobenius number 43 ; thus, for every 
77 > 43 , one can purchase exactly n McNuggets using the standard package sizes of 
6 , 9 , and 20 . As with the numerical monoid (771,772), almost every element in M has 
multiple distinct factorizations. For example, 60 can be expressed as three copies of 
20 , ten copies of 6 , or four copies of 6 plus four copies of 9 . Each such expression 
corresponds to a way to purchase exactly 60 McNuggets. 

The existence of a Frobenius number for a numerical monoid S guarantees that past 
some value, every natural number lies in S. Therefore, as a set, a numerical monoid 
S' is a scattering of “small” elements (those less than F{S)) and an infinite ray after 
F{S). What is less understood is the factorization theory of these elements, which is 
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hinted at in Example |1.4[ In general, the factorization theory of nnmerical monoids 
with 2 minimal generators is well-nnderstood, while those nnmerical monoids with 3 or 
more irredncible elements have a mnch more complicated factorization strnctnre with 
nnmerons open qnestions. 

Factorization strnctnres are often classified and qnantified using a variety of impor¬ 
tant factorization invariants, which generally measure how far from unique a monoid 
element’s factorizations are. The coming sections introduce several standard invari¬ 
ants, together with results and open questions in the context of numerical monoids. 
Section describes the strongest (but most cumbersome) invariant, the set of fac¬ 
torizations. Then, Section describes the length set and elasticity invariants, which 
measure, in different ways, how many irreducible elements appear in the factorizations 
of an element. Section gives several important results for delta sets, which measure 
the gaps in lengths that are possible between distinct factorizations of a single element. 
While length sets, elasticity, and delta sets record only information about the length 
of factorizations of an element, the catenary degree, which is described in Section 
measures how many irreducible elements different factorizations of the same element 
do not share in common. Lastly, as non-unique factorizations coincide with the pres¬ 
ence of irreducible elements that are not prime. Section describes cn-primality, which 
measures how “far from prime” a monoid element is. 


2. Sets of factorizations 

Definition 2.1. Fix a numerical monoid S = {ni,... ,nk) and an element n G S'. 
A factorization of n is an expression n = Ui of n as a sum of irreducible 

elements Ui,..., of S'. Collecting like terms, we often write factorizations of n in the 
form n = aiUi (see Remark 1.2). Write 

Z5(ii) = {(ai,..., Ofc) : n = aiUi H-h akUk] C 

for the set of factorizations of n G S'. When there is no ambiguity, we often omit the 
subscript and simply write Z{n). 


Remark 2.2. The notation for factorizations in Definition 2.1 is motivated in part by 
connections to discrete geometry. In particular, each factorization of n corresponds to 
an integer point in and the set '7.sin) forms a codimension-1 linear cross-section 
of N*’ that coincides with the set of integer points of a convex integral polytope. This 
viewpoint allows the use of techniques from lattice point enumeration and other tools 
from combinatorial commutative algebra when examining the factorization structure 
of numerical monoids 


The set of factorizations (Definition 2.3) is an example of a perfect factorization 
invariant in that it encapsulates enough information to uniquely determine the under¬ 


lying monoid (Theorem 2.4). However, such complete information comes at a cost: 
extracting information from it (or even simply writing it down) is a nontrivial task. 
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Invariants derived from the set of factorizations, snch as the length set and elasticity 
invariants (Section]^, are more manageable and easier to work with, bnt necessitate 
a loss of information. Below, 2^ denotes the power set of 


Definition 2.3. Given a nnmerical monoid S = {ni,... ,nk), let 

Z{S) = {Zs(n) ■.neS}c2^" 
denote the set of factorizations of S. 


Theorem 2.4 follows from the fact that a minimal presentation of a numerical monoid 
S (as well as its dehning congruence) can be recovered from ^(*5'). Example 2.5 ex¬ 
hibits one possible method of doing so; see [IHl SS] for more background on monoid 
congruences and minimal presentations. In fact, for a hxed n S> 0, the set Zs{n) 
alone contains enough information to recover the monoid structure of S. As such, it is 
perhaps not surprising that Z{S) is often hard to describe in complete detail. 


Theorem 2.4. For two numerical monoids S and S', the following are equivalent. 

(a) S = S', 

(b) S^ S', and 

(c) Z{S) = Z{S'). 

Proof, [ini Theorem 1.1]. □ 


Example 2.5. Fix a numerical monoid S, and suppose that the two-element sets 
{(3, 0,0), (0, 2 , 0)} and {(10, 0, 0), (0, 0, 3)} are each a subset of some element of Z{S). 
Since each is a subset of S must have three minimal generators, say S = (ni, n 2 , n^). 
Moreover, since (3, 0, 0) and (0, 2, 0) are factorizations of the same element, we must 
have 3ni = 2n2, and similarly that lOni = 3 n 3 . Since gcd(ni, n 2 , 113 ) = 1 , we conclude 


that S = (6,9,20), the McNugget monoid from Example 1.4 


Larger elements in a numerical monoid S are more likely to have larger factorization 
sets. More precisel y, if S has k minimal generators, the size of Z(n) grows on the order 
of n^~^. Theorem 2.6, which specializes the asymptotic result in [33] to numerical 


monoids, can be found in 


Theorem 2.6 ( [38] Theorem 3.9]). Given S = {ni,... ,nk) numerical, there exist 
periodic functions oq, ..., ak -2 : N —t Q whose periods divide lcm(ni,..., Uk) such that 


|Zs(n)| = 


1 


{k - l)!ni ■■■Uk 


n 


k-l 


ak- 2 {n)n H--F ai{n)n ao{n) 


for all n e S. 
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Remark 2.7. Much of the recent literature on numerical monoids has focused on 
algorithms for explicit computation. It is worth noting that although Theorem 2T 
implies |Z(?7,)| is polynomial in n, is it A^P-hard to compute this set in general, as 
doing so requires solving a system of Diophantine equations. See [22] for a survey of 
recent computational advances in this area. 


Several of the factorization invariants introduced in the coming sections, including 
those derived directly from the set of factorizations, beneht from faster computation 
methods. In particular, several of these invariants admit algorithms that do not require 
computing large sets of factorizations; see, for instance. Remarks 4.20 and 6.5 Addi¬ 


tional computational traction can also be obtained by specializing to certain families 
of monoids, such as those with a fixed number of generators or generating sets of a 


certain form; see, for instance. Theorems 3.12, 4.3 and 5.6 


Another computation technique in computing invariant values is to hrst compute 
the Apery set of S (a set determined by the complement of S'; see ffl). Although 
computing the Apery set from a list of generators is also A^P-hard in general, many 
invariants can then be computed directly from it. For instance, computing the Frobe- 
nius number of a numerical monoid is A^P-hard in general [12|, but if the Apery set 
is known, its computation is trivial jlH]. In this way, the complexity of computing an 
invariant (or several invariants) can be “transferred” to the computation of the Apery 
set |26l[35lll3|. Additionally, for some special families of numerical monoids, the Apery 
set is known to admit a closed form; see HIST]. 


3. Length sets and elasticity 

While the set of factorizations for a numerical monoid is a perfect invariant, it is 
often cumbersome to compute and encode. The next two sections focus on invariants 
obtained by passing from factorizations to their lengths. In this section, we introduce 
length sets and elasticity as useful invariants for numerical monoids. We begin with 
the dehnition of a length set. 

Definition 3.1. Fix a numerical monoid S = (ni,... ,71^), and £x n G S'. Given a 
factorization a = (oi, 02 ,.. ., a^) G Zs{n), we denote by |a| the number of irreducibles 
in the factorization a; that is, |a| = oi a?;. The set of factorization lengths of n, 

denoted L 5 (n) = {|a| : a G Zs{n)}, is called the length set of n. When there is no 
ambiguity, we often omit the subscript and simply write L(?7,). 

Example 3.2. The length set of an element is easily computed from its set of factor¬ 
izations. In the numerical monoid S' = (7,10,13), the element 20 has factorization set 
Z(20) = {(1, 0,1), (0, 2, 0)} and length set L(20) = {2}. This length set, however, is 
shared with other elements (e.g., L(14) = {2}). This example highlights one manner 
in which information is frequently lost when passing from Z(n) to L(u). 
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Although distinct elements in a numerical monoid can have identical length sets, one 
may ask whether distinct numerical monoids can have different sets of length sets. In 
what follows, 2 ^ denotes the power set of N. 


Definition 3.3. For a numerical monoid S, its set of length sets is given by 

/:(A) = {L(n) ■.neS}c2^. 

Investigations into whether sets of length sets determine finitely generated, cancella- 
tive, commutative monoids up to isomorphism are widespread in factorization theory. 
The most famous conjecture is stated below for certain block monoids (that is, monoids 
of zero-sum sequences of abelian group elements under concatenation). Since a detailed 
introduction to block monoids would take us too far afield, we refer the interested reader 
to [ 30 ] for background on block monoids and [311 [ 32 ] for recent progress on Conjec¬ 
ture 3.4 In what follows, let B{G) denote the block monoid over a group G; C{B{G)) 


is defined analogously for this class of monoids. Note that for two abelian groups G 
and G', B{G) ^ B{G') if and only if G ^ G'. 

Conjecture 3.4. Given two finite abelian groups G and G' with |G|, |G'| > 3, we have 
C{B{G)) = C{B{G')) implies B{G) ^ B{G'). 


Theorem 3.5 answers an analogous question for numerical monoids in the negative. 

Theorem 3.5. [21 Corollary 3.4] There exist numerical monoids S S' with C{S) = 
C{S'). In particular, sets of length sets do not characterize numerical monoids. 

Remark 3.6. The distinct numerical monoids with equal sets of length sets described 
in [ 2 ] are called arithmetical numerical monoids, which are numerical monoids gener¬ 
ated by an arithmetic sequence a,a + d,... ,a + kd with gcd(a, d) = 1 and 1 < fc < a. 


Theorem |3.5| is proven by computing the set of length sets of arithmetical numerical 
monoids in terms of the parameters a, d and k and showing that different values of 
these parameters yield the same set. 


The description of the set of length sets for arithmetical numerical monoids upon 
which Theorem |3.5| is based is rare and, in general, it is difficult to obtain such a descrip¬ 
tion for other families of numerical monoids. Part of the difficulty comes in computing 
every factorization length that an element in a numerical monoid can obtain. Many 
times, the most interesting factorization lengths for an element are the largest and 
smallest ones. The quotient of these two values gives a weaker factorization invariant 
called elasticity, which is defined below. 

Definition 3.7. Let S' be a numerical monoid, n E S nonzero, and L(?t,) its length set. 
Let M(n) = maxL(n) and m(u) = minL(?T,). The elasticity of n is given by 

p{n) = M(n)/m(n). 

The set of elasticities of S is given by R{S) = {pin) : n G S, n 7 ^ 0}. The elasticity of 
S is given by p(S) = sup R(S). 
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Figure 1. A plot showing the elasticities for the numerical monoid 
(20,21,45) (left) and the arithmetical numerical monoid (7,10,13,16) 
(right), produced using SAGE [19]. Here, the point at {n,q) indicates 
that p{n) = q. 


In what follows, we see that p{S) < oo for numerical monoids and that, in fact, 
p{S) = ma.x R{S). Moreover, elasticity is one of the few factorization invariants known 
to admit a closed form for any numerical monoid S. 


Theorem 3.8 ([T9l Theorem 2.1 & Corollary 2.3]). Let S = {ni,n 2 ,... ,nk) be min¬ 
imally generated by ni ,... ,nk with ni and Uk the smallest and largest generators, re¬ 
spectively. Then, p{S) = piriiUk) = Uk/ni and the set of elasticities R{S) has Uk/ni 
as its unique accumulation point. 


Example 3.9. Figure depicts the elasticity of elements in the numerical monoids 
S = (20,21,45) and S' = (7,10,13,16). The maximal values attained are 9/4 and 
16/7, respectively, as ensured by Theorem 3.8 Moreover, the monoid elements with 


maximal elasticity are precisely the multiples of the least common multiple of the hrst 
and last generators. 


More can be said about the structure of the set R{S) upon examining the asymptotic 
behavior of the maximum- and mininum-length factorization functions max 15 ( 77 ) and 
min 15 ( 77 ) for large monoid elements 77 G S'. See [6l Section 4] for more detail. 


Theorem 3.10. [H Corollary 4.5] Fix a numerical monoid S = ( 771 , 772 ,..., 77^). 
(a) For n > Uk-iUk, we have 


p{n + 77i77fc) = 


M{n) + Uk 


m{n) -|- 77i 

(b) The set R{S) is the union of a finite set and a collection of niUk monotone in¬ 
creasing sequences, each converging to rik/ni. 
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Remark 3.11. Both plots in Figure [T] demonstrate the claims in Theorem 3.10 
Roughly speaking, each plot can be viewed as a union of “wedges” of points, where 
each wedge contains one point from each of the monotone increasing sequences de¬ 
scribed in Theorem 3.10[ b). Moreover, the right-hand plot depicts the elasticies of 
an arithmetical numerical monoid, where elasticity values are often repeated in quick 
succession. It is this repetition of elasticitiy values that gives rise to Theorem |3.12 


Certainly, two numerical monoids with identical sets of length sets have equal sets 
of elasticities. Conversely, there exist pairs of numerical monoids with equal elas¬ 
ticity sets and unequal sets of length sets. For example, for S = (6,10,13,14) and 
S' = (6,11,13,14), a computation shows that R{S) = R{S'), while {4,6} G C{S) 
but {4,6} ^ R{S'); see P, Example 3.11]. Thus, in general, sets of elasticities do not 
determine sets of length sets. Theorem |3.12 demonstrates that, within the class of 
arithmetical numerical monoids, equality of sets of length sets is equivalent to equality 
of sets of elasticities. 


Theorem 3.12. [HI Theorem 1.2] For two distinct arithmetical numerical monoids 
S = {a,a + d,... ,a + kd) and S' = {a',a' + d',... ,a' + k'd'), the following are equivalent: 

(a) R{S) = R{S') 

(b) C{S) = C{S') 


Theorem 3.10 demonstrates that computing the set of elasticities is fairly tractable, 
while the set of length sets is, in general, much more difficult to describe. Therefore, 
finding classes of numerical monoids (beyond the arithmetical case) for which the set 
of elasticities is as strong an invariant as the set of length sets is of much interest. 


Problem 3.13. Characterize which distinct numerical monoids S and S' satisfy R{S) = 
R{S') and C{S) ^ L{S'). 


4. The delta set 


In this section, we consider the delta set invariant (Definition 4.1). The delta set 
of a monoid element n, like its elasticity (Definition 3.7), is derived from the set of 
factorization lengths L(n). 


Definition 4.1. Fix a numerical monoid S = (ni,... ,nk), and fix n G S'. Writing 
L(n) = {ii < ■ ■ ■ < £r}, the delta set of n is the set /S.s{n) = {li — : 2 < i < r} 

of successive differences of factorization lengths of n, and A(S') = A 5 (n). When 

there is no ambiguity, we often omit the subscript and simply write A(n). 


Example 4.2. Consider the numerical monoid S = with two relatively prime 

generators. If n G S', then there exist a,b & N such that n = ani + bn 2 - Notice that 
for any fc G Z, n = (a -f kn 2 )ni + {b — kni)n 2 is a factorization so long as a -|- kn 2 
and b — kni are non-negative. In fact, it can be shown that all factorizations of n in 
S are of this form, corresponding to successive values of k. Since the length of such a 
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factorization is a + h + k{n 2 — Ri), all successive length differences in L(n) yield n 2 — ni. 
Thus, for any n E S, either A(?7,) is empty (when its length set is a singleton), or 
A(?7,) = {n 2 — rii}. In particular, 

A(S') = {na - ni}. 

The above example is a special case of arithmetical numerical monoids, whose mini¬ 
mal generators form an arithmetic sequence (Remark |3.6| ). Such monoids have single- 
ton delta sets consisting of the step size between these generators, a fact which plays 
a crutial role in characterizing their length sets (Theorems |3.5[ ). 

Theorem 4.3. |21 Corollary 2.3] Fix positive integers a, d, and k with gcd(a, d) = 1 
and 1 < k < a. For S = {a,a + d,a + 2d,... ,a + kd) and n E S, either A(n) is empty, 
or A(n) = {d}. In particular, A(5) = {d}. 

Even slight generalizations of arithmetical numerical monoids have elusive delta sets; 
this leads us to the following open problem. 


Problem 4.4. Let a, h, d, and k he positive integers such that gcd(a, d) = 1 and k < a, 
and consider S = (a, ah + d,ah + 2d,... ,ah + kd). Compute A(S'). 


Remark 4.5. The generators of the numerical monoid in Problem 4.4 form what is 
known as a generalized arithmetic seguence. Numerical monoids generated by gener¬ 
alized arithmetic sequences and arithmetic numerical monoids are closely connected: 
when h = 1 , a generalized arithmetic sequence reduces to an arithmetic sequence; also, 
the last k generators of a generalized arithmetic sequence form an arithmetic sequence. 
These relationships, together with a precise membership criterion m Theorem 3.1], 
make working with this family of numerical monoids more tractable. 


Remark 4.6. Aside from the arithmetical setting, few natural families of numerical 
monoids are known to admit closed forms for their delta sets. In general, given a 
numerical monoid S, it can be hard to prove that d ^ A (S') for a given integer d, 
especially if min A(S') < d < maxA(S'), since it must be shown that d does not lie in 
the delta set of even a single element of S. Example |4.7 demonstrates the sublety of 
this question. 

Example 4.7. The delta set of elements of the numerical monoid S = (17,33, 53, 71) 
are depicted in Figure]^ Notice that A(S') = {2,4,6}, but 6 only appears in A(266), 
A(283), and A(300). Indeed, for n > 300, A(n) C {2,4}; see Theorem 4.16 


Remark 4.8. In contrast to numerical monoids, the delta set of any block monoid 
M = B{G) over a finite abelian group G is known to be an interval, and recent work 
has focused on determining its maximal element [301 Section 6.7]. In particular, any 
positive integer d < maxA(M) must also lie in A(M). Roughly speaking, block 
monoids have abundant irreducible elements, so if d G A (a) for some a E M, one can 
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Figlfre 2. A plot showing the delta sets of elements in the numeri¬ 


cal monoid (17,33,53,71) from Example 4.7, produced using SAGE [49] . 

Here, each point {n,d) indicates that d G A(n). 

produce an element in M with d — 1 in its delta set by simply adding the appropriate 
irreducible elements to a. This process effectively “cuts down” a length gap d by 
carefully hlling in a particular factorization length. 


The difficulty in controlling the phenomenon discussed in Remark 4A when con¬ 
sidering numerical monoids motivates Problem |4.9 , which is known as the delta set 
realization problem for numerical monoids. 

Problem 4.9. Let D C Z>o be a finite set of positive integers. Determine if there 
exists a numerical monoid S such that A(S') = D. 

Currently, the only known restrictions on A (S') for a numerical monoid S are given 


in Theorem 4.10, each of which also holds for a much larger class of monoids. 


Theorem 4.10 f [30l Proposition 1.4.4]). If M is an atomic, cancellative, commutative 
monoid with ^{M) ^ 0, then, min(A(M)) = gcd(A(M)). If, additionally, M is finitely 
generated, then A{M) is finite. 


Remark 4.11. Progress on Problem 4.9| has come from Ending specific subsets that 
can be realized as delta sets. For example, every set of the form {d, 2d ,..., td} is 
known to occur as the delta set of a numerical monoid m- Additionally, a family of 
numerical monoids whose delta sets are {d, td} for any t, d > 1 is given [22] . 


Recent investigation into Problem 4.9| has greatly benefited from evidence produced 
by computer algebra systems. That said, until very recently, there was no known 
algorithm to compute the delta set of an entire numerical monoid S. Although the delta 
set A(n) of an element n E S can be computed by computing its set of factorizations 
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Z(n), it is less straight forward to compute the entire set A (S'), as this requires taking 
the union of the delta sets of all (inhntely many) elements of S. 

Before discussing how this issue can be resolved (see Corollary 4.17), we state The¬ 
orem 4.14, which demonstrates that the maximal element maxA(S) in the delta set 


of a numerical monoid S can be obtain by computing the delta sets of a special (com¬ 
putable) class of monoid elements, called Betti elements (Dehnition 4.12). 


Definition 4.12. Let S = (ni, n 2 , • • •, R-fc) be a numerical monoid and n G S'. Consider 
the graph Gn with vertex set given by its set of factorizations Z(n) and an edge connect¬ 
ing a, b G Z(n) if a and b have disjoint support as vectors. That is, a = (oi, ..., a^) 
and b = {bi,... ,bk) are adjacent in Gn if for all i, ai and bi are never both non-zero. 
An element n G S' is called a Betti element if Gn is disconnected. The set of Betti 
elements of S is denoted by Betti(S'). 


Remark 4.13. Every numerical monoid contains only hnitely many Betti elements. 
When S' = ( 17 . 1 , 112 ) is a numerical monoid with two irreducible elements, nin 2 is the 
unique Betti element. In the case when S' = (ni, n 2 , nf) has three irreducible elements, 
there are at most three Betti elements, each of which is an integer multiple of one of 
the generators of S' m Section 9.3]. 


Theorem 4.14 f [THl Theorem 2.5]). Let S be a numerical monoid. The largest element 
0 /A(S') is achieved at a Betti element of S. In other words, 

max A(S') = max{d : d G A{b) for b G Betti(S')}. 


Betti elements are discussed in more depth in Section as they are useful in finding 
maximal and minimal catenary degrees achieved in a numerical monoid (see Defini¬ 


tions 5.2 and 5.4). 


Example 4.15. As stated in Remark 4.13, for the numerical monoid S' = (ni, 772 ) with 


two irreducible elements, 77i?72 is the unique Betti element. In fact, it has precisely two 
factorizations, namely {n 2 , 0) and (0, ni). Since their difference in length equals n 2 —ni, 


we have max A(S') = 112 — rii. Of course, by Theorem 4.3, A(S') = {112 — rii}. 


Theorem 4.16 states that delta sets for numerical monoids are eventually periodic 
with period dividing lcm(77i, rifc). The value of Ns given in [25] is explicitly computed. 


albeit cumbersome to state. Most impressive about Theorem |4.16| is that it provides a 
concrete upper bound for the beginning of the periodic behavior, and thus a hnite set 
of monoid elements the union of whose delta sets equals the delta set of the entire nu¬ 


merical monoid. As a consequence. Theorem 4.16 also yields the hrst known algorithm 


to compute the delta set of any numerical monoid; see Corollary 4.17 


Theorem 4.16 f [25l Theorem 1]). For S = {ni,... ,nk) numerical, there exists a 
number Ns such that for any n ^ S with n > Ns, A(n) = A(?7 -|- lcm(?7i, Ufc)). 
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Corollary 4.17. Given a numerical monoid S = (ni,n 2 ,... ,nk), we have 

A(S)= U A(n), 

nesn[o,A] 

where N = Ns + lcm(ni,nfc). 



for all n > 319. 


Remark 4.19. While Corollary 4.17 provides the delta set of the entire numerical 
monoid by taking the union of only hnite many element-wise delta sets, the number of 
elements n for which A(n) needs to be computed becomes prohibitively large, especially 
when larger generators are used. Aside from recent improvements in the special case 
when S has 3 minimal generators izu, few refinements are known for the set of elements 
whose delta sets must be computed. Alternatively, we can utilize a dynamic algorithm 


for computing length sets to dynamically compute A(n) and A(S'); see Remark 4.20 


Remark 4.20. In the numerical monoid S = (6, 9,20), it is possible to compute L(60) 
directly from L(40), L(51), and L(54), without computing Z(60). The key observation 
is that one can produce a factorization of 60 from a factorization of 40, 51 or 54 by 
adding a single irreducible element, and this opereration always increases factorization 
length by exactly one. Consider the following computation: 

L(60) = (L(40) + 1) U (L(51) + 1) U (L(54) + 1) 

= {3} U {7,8,9} U {7,8,9,10} 

= {3,7,8,9,10}. 

This provides an inductive (dynamic) method to compute all the length sets necessary 
to employ Corollary |4.17| without having to compute any factorization sets, and the 
resulting algorithm is signihcantly faster than one dependent on factorizations. Indeed, 
the set Z(?T,) of factorizations of n grow large very quickly [33], whereas |L(?t,)| grows 
linearly in n [6l Theorem 4.3]. For a more detailed analysis, see [71 Section 3]. 


5. The catenary degree 

Many of the factorization invariants discussed above (e.g., length sets, delta sets, 
and elasticity) measure only information about the lengths of factorizations of an ele¬ 
ment. While helpful in measuring the non-uniqueness of factorizations of an element, 
these invariants lose certain information when passing from the set of factorizations 
to the set of lengths (Theorem 3.5). The catenary degree, which we describe in this 
section, is derived directly from factorizations and is not uniquely determined by sets 
of factorization lengths. 
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Definition 5.1. Fix a numerical monoid S = {rii,... ,nk), and fix an element n E S. 
For a, b G Z(n), the greatest common divisor of a and b is given by 

gcd(a, b) = (min(ai, 6 i),..., min(ar, br)) G 

and the distance between a and b (or the weight of (a, b)) is given by 

(i(a, b) = max(|a — gcd(a, b)|, |b — gcd(a, b)|). 

The function d : Z(n) x Z(n) —?■ N measures how many irreducible elements two 
factorizations do not share in common. As implied by the notation, d{—, —) is a metric 
on the set Z(n) m Proposition 1.2.5.3]. In fact, it is common to represent Z(n) as 
the vertices of a complete graph, where the edge between factorizations a and b is 
decorated with d{a, b). 

Below, we use the metric d to define the catenary degree, which measures the “overall” 
distance between distinct factorizations of a given monoid element. 

Definition 5.2. Given a, b G Z(n) and > 1, an N-chain from a to b is a sequence 
ai,...,ar G Z(n) of factorizations of n such that (i) ai = a, (ii) a^ = b, and (in) 
(i(aj_i,aj) < N for all i < r. The catenary degree of n, denoted c{n), is the smallest 
N eN such that there exists an A^-chain between any two factorizations of n. 


With respect to the decorated complete graph representing Z(n), the catenary degree 
c(n) can be computed by repeatedly removing edges of maximal distance until the graph 
becomes disconnected. The distance decorating the last removed edge is equal to c(n). 
Example |5.3| describes such a computation. 


Example 5.3. Consider the numerical monoid S = (11, 36, 39). The left-hand picture 
in Figure [^depicts the factorizations of 450 G S along with all pairwise distances. There 
exists a 16-chain between any two factorizations of 450; one such 16-chain between 
(6,2,8) and (24,3,2) is depicted with bold red edges. Since every 16-chain between 
these factorizations contains the edge labeled 16, we have c(450) = 16. 

This can also be computed in a different way. In the right-hand picture of Figure 
only distances of at most 16 are depicted, and the resulting graph is connected. Remov¬ 
ing the edge labeled 16 yields a disconnected graph, so we again conclude c(450) = 16. 


While the catenary degree of an individual element contains information about its 
different factorizations, the range of catenary degrees achieved by all elements in a 
monoid has also been studied. Definition 5A gives two ways to measure how the 
catenary degree varies throughout a numerical monoid. 


Definition 5.4. The set of catenary degrees of S is the set C(S') = {c(m) : m E S'}, 
and the catenary degree of S is the maximal catenary degree: c(S') = maxC(S'). 


By far, the more well-understood of these two invariants is c{S), which has been 
computed for several classes of numerical monoids. Theorem 5.5 gives a summary of 
some known results (additionally, see Theorem 5.8). 
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( 0 . 6 . 6 ) ( 27 . 1 . 3 ) 



Figure 3. Computing the catenary degree of 450 G S' = (11,36,39), 
as in Example 5.3 Each vertex is labeled with an element of Z(450), 


and each edge is labeled with the distance between the factorizations at 
either end. The diagram on the left depicts all edges, and the diagram 
on the right includes only those edges labeled at most c(450) = 16. Both 
graphics were created using the computer algebra system SAGE |19] . 


Theorem 5.5 f JSTl 137] ). Fix a numerical monoid S. 

(a) If S = {ni,n 2 ), then c(S') = n 2 . 

(b) If S = {a,ah + d,... ,ah + kd) for positive integers a, h, fc, d with gcd(a, d) = 1 and 
k < a, then c(S') = [|] h + d. 

(c) If S = {{]Tj=i^j){Il’j=i+i(^j) : 0 < i < A;), where 2 < ai < bi and gcd{ai,bj) = 1 
for all i ^ j < k, then c(S') = max{6 i,..., bk}- 

While a closed form for c(n) for every element n is a numerical monoid is rare, such 
a formula has been found for numerical monoids generated by an arithmetic sequence. 


Theorem 5.6 ([121 Theorem 3.1]). Suppose S = (a, a + d, 
integers a, fc, d with gcd(a, d) = 1 and k < a. Then 


, a + kd) for positive 


{ 0 

|Z(n)| = l 

\ ^ 

L(n) = 1 and 

1 c(^) 

L(n) > 1 


for each nonzero n E S, where c{S) = [|] + d. In particular, C(S') = {0,2,c(S')}. 

While it remains a lofty goal to give a closed form for the catenary degree of individual 
elements in the general setting, the set of catenary degrees may prove to be more 
approachable. Problem 5.7 is, in many ways, an analog of Problem 4.9 for the set of 
catenary degrees. 


Problem 5.7. Determine which finite sets egual C{S) for some numerical monoid S. 
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Figure 4. A SAGE plot [19] showing the catenary degrees for the nu¬ 


merical monoid S = (11, 25, 29) discussed in Example 5.9 


Many of the arguments that compute the maximal catenary degree c(S') for a finitely 
generated monoid S focus on the Betti elements of S (Definition |4.12[ ). The factoriza¬ 
tions of Betti elements contain enough information about Z(S') to give sharp bounds 
on the catenary degrees occuring in C(S') (Theorem 5.8), and can be readily computed 
using the GAP package numericalsgps 


Theorem 5.8 f [T71 HT] b For any finitely generated monoid S, we have 

c{S) = maxC(S') = max{c(6) : b G Betti(S')} and 
minC(S') = min{c(&) ; b G Betti(S')}. 

Example 5.9. Let S = (11,25, 29). The catenary degrees of S are plotted in Figure]^ 
The only Betti elements of S are 58, 150, and 154, which have catenary degrees 4, 12, 
and 14, respectively. The Betti elements can be seen in the plot as the hrst occurence 
of each of these three values. Notice that c(175) = 11 is distinct from each of these 
values. By Theorem |5.8| , every element of S with at least two distinct factorizations 
has catenary degree at least 4 and at most 14. 


As with the delta set, understanding the long-term behavior of the catenary de¬ 
gree (when viewed as a function c : S' —>■ N) is of particular interest. Theorem 5.10 


shows that the catenary degree in numerical monoids is also eventually periodic, a 
phenomenon that can be seen in Figure]^ 


Theorem 5.10 ([I2l Theorem 3.1]). Fix a numerical monoid S = (ni,... ,nk)- The 
catenary degree function c : S ^ N is eventually periodic, and moreover, its period 
divides lcm(ni,... ,nk). 


Remark 5.11. Among the eventual periodicity results in this paper (Theorems 3.10 
4.16, 5.10, and 6.9), Theorem 5.10 provides the most room for improvement. There 


is currently no known bound on the start of this periodic behavior, as its proof in 
relies on the fact that a nonincreasing sequence of positive integers is eventually 
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constant. Additionally, the period is only known to divide the least common mnltiple 
of the monoid generators, but in practice this bound is far from sharp. It remains an 
interesting problem to refine understanding of these properties. 


Examination of the delta set and the cn-primality invariant (intoduced in Section 

and|6.5 An answer 


both benefit from the use of dynamic algorithms; see Remarks 4.20 


to Problem 5.12 is likely within reach, and would aide in rehning Theorem |5.10 


Problem 5.12. Find a dynamic algorithm to compute the catenary degree of numerical 
monoid elements. 


6. W-PRIMALITY 


In cancellative, commutative monoids, non-unique factorizations arise from irre¬ 
ducible elements that are not prime. The rich factorization theory in numerical monoids 
stems in part from the fact that no element is prime. In this section, we develop the 
notion of ca-primality, which measures how “far from prime” an element is in a natural 
way (see PSO] for a more thorough introduction). We begin with a definition. 


Definition 6.1. Let S = {ni,n 2 ,... ,nk) be a numerical monoid. For n ^ S, define the 
u-primality of n to be u:{n) = m if m is the smallest positive integer with the property 
that whenever a G satisfies — n E S with |a| > m, there exists a b G 

with bi < Qi for each i such that Yli=i — n E S and |b| < m. 


In a multiplicative monoid M, an element n E M is prime if whenever n \ ab for a,b E 
M, then either n | a or n | 6. Definition 6.1 is a generalization of primeness, written 
additively in the context of a numerical monoid. In fact, in a general cancellative, 
commutative monoid, Lj{n) = 1 if and only if n is prime. 

In practice, the ca-primality of an element n of a numerical monoid S = (rii,..., n^) is 
computed by Ending those a G that are maximal amo ng those with the property that 
Yl!i=i cani — n lies in S (Definition |6.2|). Proposition 


to the ca-value of n. 


6.3 


states their precise relationship 


Definition 6.2. We say that a G is a bullet for n if (i) (X]r=i — n E S and 
(ii) (Ell aini)—nj — n ^ S whenever Oj > 0. A bullet a for n is maximal if |a| = u:{n). 
The set of bullets for an element n is denoted by bul(n). 


In general, an element n has several bullets with different lengths, 
the largest such length coincides with the u:{n 


Proposition 6.3 


As stated in 
justifying the above 


definition for a maximal bullet. 


Proposition 6.3 f [4ni Proposition 2.10]). Given any numerical monoid S, 

u{n) = max{|a| : a G bul(n)} 


for each element n E S. 
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Example 6.4. Let S = (6,9,20) denote the numerical monoid from Example 1.4 
Consider the following sets: 

bul(40) = {(0,0, 2), (4,4,0), (7, 2,0), (10, 0,0), (1,6,0), (0,8,0)}, 
bul(51) = {(0,7,0), (10, 0,0), (4, 3,0), (1,5,0), (0,0, 3), (7,1,0)}, 
bul(54) = {(9, 0,0), (6, 2,0), (0,6,0), (3,4,0), (0,0, 3)}, 
bul(60) = {(4,4,0), (7, 2,0), (10, 0,0), (1,6,0), (0,8,0), (0,0, 3)}. 

We see that (10, 0, 0) G bul(40) since 10 ■ 6 — 40 = 20 G S', but 9 ■ 6 — 40 = 14 ^ S. 
Similarly, (7,2, 0) G bul(60) since 7-6 + 2- 9 — 60 = 0 gS, but omitting a single 6 or 
9 produces a difference lying outside of S. Once the above sets have been determined. 


Proposition 6.3 implies that a;(54) = 9 and a;(40) = a;(51) = a;(60) = 10. 


Remark 6.5. Several recent works examine methods of computing ca-primality in nu¬ 
merical monoids 0 El [ini EH- At the time of writing, the fastest known algorithm 
(which is implemented in the GAP package numericalsgps [23]) utilizes dynamic pro¬ 
gramming to compute bullets [7], similar to the dynamic algorithm used to compute 
the delta set (Remark |4.20[ ). Example |6.6| demonstrates the idea behind this algorithm. 

Example 6.6. Resuming notation from Example |6.4| for each bullet a G bul(54), 
either a G bul(60) or a -|- ei G bul(60). Notice that it is impossible for both of these 
to lie in bul(60) by Definition |6.2[ ii). Similarly, bullets in bul(51) and bul(40) give 
rise to bullets in bul(60). Moreover, each bullet for 60 is the “image” of a bullet for 
54, 51 or 40 in this way. Much like the recurrence relation for sets of factorizations 
(Remark [4.20 ), the set of bullets of an element n G S' is determined by the “images” 
of the bullets of n — 6, n — 9 and n — 20; see [7| for more detail. 

In general, the existence of non-unique factorizations coincides with the existence 
of non-prime irreducible elements. As such, the cn-primality invariant was originally 
constructed to measure how far irreducible elements were from being prime (see, for 
example, BIB). Although some progress has been made in the setting of numerical 
monoids 13 HI [21], Problem |6.7| remains largely open. 


Problem 6.7. Determine the u-values of numerical monoid generators. 


Remark 6.8. Both the catenary degree (Definition 5.2) and a variation called the 


tame degree share a surprising connection to ca-primality. Given n G S', the tame degree 
t(n) equals the minimum value m such that every factorization of n with at least one 
zero entry is distance at most m (in the sense of Definition |5.1 ) from a fully supported 
factorization of n, or t(n) = 0 if n has no fully supported factorizations mm- As with 
the catenary degree, the tame degree of S' is defined as t(S') = sup{t(n) : n G S'} and is 
known to be finite. It is not hard to show that c(S') < t(S'), but what is perhaps more 
surprising is that c(S') < Ci;(S') < t(S'), where Ci;(S') = max{ci;(ni),..., a;(nfc)} is the 
maximum cj-value obtained at an irreducible HD]. This connection has been explored 
in some of the recent literature, in part as a method of bounding a;(S') (TUI ESI EH]- 
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Figure 5. A plot showing the w-function for the numerical monoids 
(4,13,19) (left) and (6,9,20) (right), produced using SAGE |49] . 


Recent investigations focusing on the w-values of all non-unit elements in monoids 
have uncovered interesting asymptotic behavior. In the setting of numerical monoids, 
cj-primality grows in a periodic, linear fashion when viewed as a function u : S ^ N. 


Theorem 6.9 ([TJ [2ll [39]). Let S = (ni, ^ 2 ,..., rzfc) be a numerical monoid. The 
u-primality function is eventually quasilinear. More specifically, for Nq = we 

have uj{n) = -^n -|- ao(n) for every n > Nq, where oq : N —)• Q is ni-periodic. 


Theorem 6^ implies that the graph of the function a; : S' —)■ N eventually has the 
form of Hi lines with common slope 1/ni. Figure [^demonstrates this phenomenon. 


Remark 6.10. We resume notation from Theorem |6.9[ In contrast to Theorem 4.16 
the period of the function oq is known to be exactly ni [7|. Additionally, the value 
of Nq, though not tight, is in general a relatively good bound on the start of the 
quasilinear behavior of us [3 Remark 6.8]. Each of these improvements on the original 

in 


statement of Theorem 6.9 


are heavily motivated by data produced by the 


dynamic algorithm described in Remark 6.5 
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